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We discuss recent observations of high energy cosmic ray positrons and electrons in the context 
of hadronic interactions in supernova remnants, the suspected accelerators of galactic cosmic rays. 
Diffusive shock acceleration can harden the energy spectrum of secondary positrons relative to that 
of the primary protons and electrons and thus explain the rise in the positron fraction observed by 
PAMELA above 10 GeV. We normalize the hadronic interaction rate by holding pion decay to be 
responsible for the gamma-rays detected by HESS from some SNRs. By simulating the spatial and 
temporal distribution of SNRs in the Galaxy according to their known statistics, we are able to then 
fit the electron (plus positron) energy spectrum measured by Fermi. It appears that IceCube has 
good prospects for detecting the hadronic neutrino fluxes expected from nearby SNRs. 
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I. INTRODUCTION 

Recently, the PAMELA collaboration has published 
data on the positron fraction, Jg+/(Jg- + Je+), in galac- 
tic cosmic rays (OCR) which is seen to increase between 
~ 10 and 100 GeV [l], in contrast to the prediction of 
the standard GCR propagation model [2 . The model 
assumes that the positrons are secondaries created via 
interactions of GCR protons and nuclei with interstel- 
lar matter, hence their spectrum should be softer than 
that of primary electrons and the positron fraction should 
thus decrease with energy [3 . The combined differential 
flux of GCR electrons and positrons, (Jg- + Je+)i has 
also been measured with the Fermi Large Area Telescope 
(LAT) and is approximately fitted by a power-law 
in energy up to ~ 1 TeV [4J. Measurements by the HESS 
collaboration [5] [6] show significant steepening of the 
spectrum beyond ~ 1 TeV, while agreeing well with the 
Fermi data at lower energies. Although the Fermi LAT 
data does not confirm the sharp feature claimed earlier 
by the ATIC collaboration [7] , there does appear to be a 
small excess flux above ~ 100 GeV in comparison with 
standard GCR propagation models [UIH]. 

Both these experimental findings have generated a lot 
of interest because they may be an indirect signature of 
dark matter particles. Annihilation or decay of galactic 
dark matter can produce electrons and positrons with 
a spectrum considerably harder than that of primary 
electrons. Besides the fine-tuning challenges such mod- 
els face [HI |9] , other cosmic ray data provide important 
constraints. The antiproton-to-proton ratio observed by 
PAMELA (To] is in good agreement with the predic- 
tion of secondary production by GCRs and thus rules 
out most dark matter annihilation/decay models which 
have hadronic final states. Even purely leptonic annihi- 
lation channels are strongly constrained by the Galactic 
synchrotron radio background [11] and by the Galactic 
gamma-ray background [T? . In fact pulsars may produce 
a hard spectrum of electron-positron pairs in the right en- 
ergy range to explain both the positron flux anomaly and 
the observed electron flux [131 HI] • 



It has long been believed ^ that GCRs are generated 
by diffusive shock acceleration (DSA) [161 HZ] iii super- 
nova remnants (SNRs) p!8l. Hadronic interactions of the 
accelerated protons will create tt^ (and tt^) which then 
decay to yield secondary and neutrinos (and 7-rays). 
It has been suggested that the acceleration of the sec- 
ondary positrons in a nearby SNR shock wave may be 
responsible for the PAMELA anomaly [19 . The fraction 
of secondary e+ which are accelerated increases with en- 
ergy, so their final spectrum is harder than the injected 
spectrum. This effect had been noted earlier as a gen- 
eral expectation for the secondary-to-primary ratio in the 
presence of stochastic Fermi acceleration [20l [21] . A sim- 
ilar effect is then predicted at higher energies for both an- 
tiprotons [22, 23 and for secondary nuclei such as boron 
[24] . These predictions will be tested soon with data from 
PAMELA and the forthcoming AMS-02 mission [25j . 

It is interesting to ask whether this model can account 
also for the absolute fluxes of e~ (and e+) in GCR at 
energies > 100 GeV [19]. This is rather sensitive to the 
assumed spatial distribution of the sources so in this pa- 
per we consider a realistic distribution of SNRs based on 
astronomical data (Sec. [ll|. Previously the flux of sec- 
ondary e~ and e+ in the sources has been normalized 
with respect to the primary electrons in an ad hoc fashion 
[T9] . Instead, we exploit the hadronic origin of these sec- 
ondaries and normalize using the 7-ray fluxes (assumed 
to be from tt^ decay) detected from known SNRs by 
HESS. We can thus fix the only free model parameter 
by fitting the total e~ + e+ flux to Fermi LAT and HESS 
data (Sec. HI). The e+ fraction is then predicted up to 
TeV energies and provides a good match to PAMELA 
data (Sec. |IV[ ). Having constrained the distribution of 
the closest SNRs via the measured e~ and e+ spectra, we 
present an example of a likely source distribution in order 
to illustrate that there are good prospects for IceCube to 
detect neutrinos from nearby SNR. A consistent picture 
thus emerges for all presently available data in the frame- 
work of the standard DSA/SNR origin model of GCR. 
However there remain some open issues and grounds for 
concern which we discuss at the end (Sec. [V]). 
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FIG. 1: Distance-time diagram for nearby SNRs (after Ref. [30]). Left: The open circles mark supernova events and the 
world-hnes of the discovered remnants are indicated. The thick yellow line is our past light-cone; all events lying on it, e.g., the 
SNR world-lines touching it, can be observed presently. The blue, purple and red shadings (top to bottom) show the relative 
contribution of sources to the diffuse e~ and flux observed at Earth at 10, 100 and 1000 GeV, respectively. The open red 
circle is an example of a hypothetical supernova whose remnant is too old to be visible any longer but which might still be 
contributing to the diffuse e~ and flux. Right: A distance-time diagram for hypothetical nearby SNRs. The open black 
circles are an example of a possible "history" of supernovae (the wordlines have been suppressed) as simulated by our Monte 
Carlo calculation (see Section 



II. DIFFUSION MODEL AND SOURCE 
DISTRIBUTION 

The diffusive transport of high energy e~ and e+ in 
the Galaxy is governed by the equation [15j, 

^ = V(Z)GCRVn±) + A(6„^)+Q^, (1) 

where n±dE = n±{r,t^ E)dE denotes the particle den- 
sity of e+ and e~ with energy in [E^E + dE]. The 
spatial diffusion coefficient is assumed not to depend 
on the position in the Galaxy but only on energy: 
^gcr(^) = DqE^. The energy loss rate of GCR 
and e+ through synchrotron radiation in Galactic mag- 
netic fields and inverse Compton scattering on the CMB 
and interstellar radiation backgrounds is parametrized as 
b{E) = boE'^. Finally, Q± denotes the injection of elec- 
trons and positrons from both (possible) primary and sec- 
ondary sources. In the majority of previous calculations, 
in particular the GALPROP code [2] in its conventional 
setup, the distribution of sources is assumed to be con- 
tinuous. However at energies when the diffusion length i 
becomes smaller than the distance to the closest source 
(similar to the average distance between sources for a 
homogeneous distribution) the fact that the sources are 
discrete should become important. This effect on GCR 
electrons was first pointed out in Ref. ^ and has later 
been considered in more detail [27l [28l |29] . It has also 
been studied using an extended version of GALPROP 
|2Qj. 



Assuming that all the electrons and positrons are re- 
leased instantly at the end of the Sedov- Taylor phase of 
expansion when the SNR becomes radiative, the flux of 
N sources at distances and times ti is given by the 
sum over the corresponding Green's functions (see Ap- 
pendix |A| for the particle density (times the usual "flux 
factor", c/47r, for an isotropic population of relativist ic 
particles): 

N 

^iv(^) = ^^Gdisk(^,r„t,)- (2) 

The spatial distribution and temporal history of GCR 
sources, {ri,^^}, is not known a priori but can be mod- 
elled for an assumed source class e.g., SNRs. At low 
energies, the GCR diffusion length is long enough such 
that the approximation of a continuous source density is 
acceptable. For energies above some hundreds of GeV, 
however, the fluctuations introduced by the discreteness 
of the sources cannot be neglected any longer. 

Some authors [28, 31 have assumed a continuous dis- 
tribution of sources for distances beyond a few hundred 
parsecs, supplemented by a set of known SNRs for smaller 
distances. This approach is however biased by the choice 
of young, nearby sources which have been detected in ra- 
dio and/or X-rays. Older sources may not be visible in 
photons any longer but still be contributing to the GCR 
electron flux, see Fig. [l] We note that the effect of this 
incomplete assumed source distribution is a dip in the 
electron flux seen in both analyses ^Hl [31] , although at 
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FIG. 2: Left: The assumed distribution of SNRs in the Galaxy; the cross denotes the position of the Sun in between two spiral 
arms. Right: The probability density for the distance of a SNR from the Sun. 



different energies because of the different diffusion model 
parameters chosen. 

Determining the complete distribution of sources in our 
vicinity {i.e. up to a few kpc) from observations seems 
challenging. However it turns out that we do not need 
to know the exact distribution in order to calculate the 
e+ flux but require only a limited amount of information, 
most of which is encoded already in the dominant e~ flux. 
By including the recent measurements by Fermi LAT [4] 



and HESS [5] |6] of the total e' 



flux in the energy 



region of interest, we have sufficient information at hand 
to make a prediction for the positron fraction under the 
assumption that the additional positrons originate in the 
same sources. 

We perform a Monte Carlo calculation by considering 
a large number of random distributions of sources drawn 
from a probability density function that reflects our as- 
tronomical knowledge about the distribution of SNRs in 
the Galaxy. The better the flux of e~ and e+ from such a 
"history" of sources reproduces the measured fluxes, the 
closer is the underlying distribution of sources likely to be 
to the actual one. Of course all SNRs are not the same, 
however variations of the source parameters would only 
introduce additional fluctuations into the fluxes without 
altering their average. We can choose the best "fit" to 
the data and thus determine the e+ flux. 

The smoothed radial distribution of SNRs in the 
Galaxy is well modelled by [32) : 



/(r) = A sin 



To 



+ ^0 e 



-/3r 



(3) 



where A = 1.96 kpc"^ ro = 17.2 kpc, 6>o = 0.08 and 
P = 0.13. To obtain a realistic probability density for 
the distance between the Earth and a SNR we have to 
also take into account the spiral structure of the Galaxy. 



We adopt a logarithmic spiral with four arms of pitch 
angle 12.6° and a central bar of 6 kpc length inclined by 
30° with respect to the direction Sun - Galactic centre 
[33 . The density of SNRs is modelled by a Gaussian 
with 500 pc dispersion for each arm [27]. The resulting 
distribution ^(r, ^) (see left panel of Fig. [2| has been 
normalized with respect to azimuth in such a way that 
the above radial distribution (|3| is recovered. To obtain 
the probability density for the distances we transform to 
the coordinates {r'^cj)') centered on the Sun. As the e~ 
and e+ fluxes are assumed to be isotropic, we can average 
over the polar angle 0', such that the probability density 
fr' depends only on the distance r' to the source. 



fr'{r') = ^ d^'r'g{r{r', ^'), ^r', ^' j). 



(4) 



This function is shown in the right panel of Fig. [2] 

We assume that the sources are uniformly distributed 
in time, i.e. their probability density ft{t) is 



ft{t) 



otherwise, 



(5) 



with tmax standing for the earliest time considered, which 
is related to the minimum energy for which our calcula- 
tion is valid through: 

-1 



(6) 



The total number N of sources that are needed in 
the Monte Carlo to reproduce the (observed) number 
JV 300 of SNRs active in the galaxy at any given 
time depends on the average lifetime of a SNR, tsnr, 



which is suggested to be 
300 



10"^ yr [TH;, hence 



TV = 3 X 10^ 



108 yr / \l0^yr 



TSNR 



(7) 
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Source 


Other name(s) 


F 


J° 


- 10-^^ 


-E'max 


d 


Q° ± 10^^ 


Ref. 








[(cm2 


sTeV)-i] 


[TeV] 


[kpc] 


[(sTeV)-i] 




HESS J0852-463 


RX J0852.0-4622 (Vela Junior) 


2.1 ±0.1 


21 


±2 


> 10 


0.2 


0.10 


m 


HESS J1442-624 


RCW 86, SN 185 (?) 


2.54 ±0.12 


3.72 


±0.50 


> 20 


1 


0.46 


m 


HESS J1713-381 


CTB 37B, G348.7+0.3 


2.65 ±0.19 


0.65 


±0.11 




7 


3.812 


m 


HESS J1713-397 


RX J1713.7-3946, G347.3-0.5 


2.04 ±0.04 


21.3 


±0.5 


17.9 ± 3.3 


1 


2.55 


[371 [38] 


HESS J1714-385 


CTB 37A 


2.30 ±0.13 


0.87 


±0.1 


> 12 


11.3 


13.3 


m 


HESS J1731-347 


G 353.6-07 


2.26 ±0.10 


6.1 


±0.8 


> 80 


3.2 


7.48 




HESS J1801-233" 


W 28, GRO J1801-2320 


2.66 ±0.27 


0.75 


±0.11 


Z 4 


2 


0.359 


m 


HESS J1804-216^ 


W 30, G8.7-0.1 


2.72 ±0.06 




5.74 


> 10 


6 


24.73 


m 


HESS J1834-087 


W 41, G23.3-0.3 


2.45 ±0.16 




2.63 


>3 


5 


7.87 


m 


MAGIC J0616+225 


IC 443 


3.1 ±0.3 




0.58 


Z 1 


1.5 


0.156 


m 


Cassiopeia A 




2.4 ±0.2 


1.0 


±0.1 


> 40 


3.4 


1.38 




J0632+057 


Monoceros 


2.53 ±0.26 


0.91 


±0.17 


N/A 


1.6 


0.279 


m 


Mean 


- 2.5 






>20 




- 5.2 




Mean, excluding sources with F > 2.8 


- 2.4 






>20 




- 5.7 




Mean, excluding sources with F > 2.6 


- 2.3 






>20 




- 4.2 





^We assume that W 28 powers only the emission from J1801 
^W30 is taken to be the origin of the VHE emission [43) . 
^Cas A was first detected by HEGRA 47 . 



-233 (and not the nearby J1800-240 A, B and C). 



TABLE L Summary of spectral parameters for SNRs detected in 7-rays from a power-law fit to the spectrum, = 
Jj{E /TeY)~^ , with an exponential cut-off at ^max in the case of HESS J1713. 7-397. The errors shown are statistical only — 
the systematic error is conservatively estimated to be 20% on the flux and ±0.2 on the spectral index F. Also shown is the 
estimated distance d and the injection rate derived from Eq. (24). 



III. FITTING THE e+ ± e' SPECTRA 



A. Primary electrons 



A schematic description of the present framework is 
shown in Fig. |3] Cosmic rays are shock accelerated in 
SNRs and then diffuse through the Galaxy to Earth un- 
dergoing collisions with interstellar matter en route and 
creating secondary e+. As discussed, the ratio of the sec- 
ondary e+ to the primary e~ from the sources should 
decrease with energy, in contrast to the behaviour seen 
by PAMELA. We follow Ref. [19] in explaining this by in- 
voking a new component of e+ which is produced through 
cosmic ray interactions in the SNRs, and then shock ac- 
celerated, thus yielding a harder spectrum than that of 
their primaries. We discuss these components in turn 
below and calculate their relative contributions by nor- 
malising to the 7-ray flux from the SNRs, which pro- 
vides an independent measure of the hadronic interac- 
tions therein. 



Acceleration in SNR 



Propagation in Galaxy 



accelerated 
secondary 
secondary (i/, 7) 




new 
component 



conventional 
component 



FIG. 3: Schematic description of contributions to the galactic 
cosmic rays observed at Earth in the present framework. 



The radio and X-ray emission observed from SNRs is 
interpreted as synchrotron radiation of electrons acceler- 
ated up to energies of O(IOO) TeV [18]. The spectrum of 
this radiation then determines the spectrum of the un- 
derlying relativistic electrons. Moreover the theory of dif- 
fusive shock acceleration [161 El predicts similar spectra 
for the accelerated protons and nuclei as for the electrons. 
If the 7-ray emission observed by HESS from a number 
of identifled SNRs is assumed to be of hadronic origin, 
we can use the measured spectra to constrain both the 
relativistic proton and electron population. 

Table [T| shows a compilation of 7-ray sources observed 
by HESS that have been identifled as SNRs. We have in- 
cluded all identifled shell-type SNRs and strong SNR can- 
didates in the HESS source catalogue [49], and also added 
the SNRs IC 443, Cassiopeia A and Monoceros. Actually 
it is not clear that the acceleration of secondaries does 
occur in all the SNRs considered, especially when the 7- 
ray emission is associated with a neighbouring molecular 
cloud rather than coming from the vicinity of the shock 
wave. In fact the 7-rays could equally well be due to 
inverse-Compton scattering by the relativistic electrons 
responsible for the observed synchroton radio and X-ray 
emission. Therefore, we have considered three possibil- 
ities — including all sources implies a mean power-law 
spectral index for the protons of (F) = 2.5, while exclud- 
ing steep spectrum sources with F > 2.8 gives (F) = 2.3 
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and excluding sources with F > 2.6 yields (F) = 2.4. In 
the following we adopt the central value, F = 2.4, for the 
electron population too, unless stated otherwise. This 
requires a compression factor of r < 3.3 in contrast to 
the value of r = 4 expected for a strong shock, so there 
is clearly some tension between the DSA theory and ob- 
servations. This can possibly be resolved if we consider 
only a subset of the SNRs in Table 1 to be hadronic ac- 
celerators, or if the 7-ray spectrum is steepened e.g. by 
the onset of an exponential cutoff in the electron spec- 
trum. Our model assumptions are intimately connected 
to the production of neutrinos, the detection of which 
will therefore provide an independent test as we discuss 
later. In this work we adopt a cut-off of ^cut — 20 TeV 
which is consistent with DSA theory [18]. The source 
spectrum of primary electrons is then: 



f— y 

VGeVy 



(8) 



The normalisation R^^_ is determined by fitting the elec- 
tron flux at Earth resulting from our Monte Carlo com- 
putation to the preliminary measurement by PAMELA 
at 10 GeV [50^; the secondary fluxes can be neglected for 
this normalisation. We find = 1.8 x 10^^ GeV"^ for 

e 

F = 2.4 which corresponds to a total injection energy of 



/ dEE R^- {E)c^7x 10^^ erg. 



(9) 



This compares well to the value of 9.2 x 10^^ erg said to 
be required to power the GCR electrons [18 . 

Solar modulation which is important below ~ 10 GeV, 
has been accounted for using the force field approach 
[ST, with a charge-independent potential of (j) = 600 MV. 
However, our simple model ignores convection and 
(re) acceleration in the interstellar medium which become 
important below ^ 5 GeV, hence the electron flux cannot 
be predicted at lower energies. The primary e~ fluxes as 
measured on Earth for 30 different source configurations 
are shown in the top panel of Fig. |5j With an injection 
power-law index F :^ 2.4 ±0.1 as required for consistency 
with the 7-ray data, there clearly is a deficit at high en- 
ergies compared to the e+ -h e~ flux measured by Fermi 
LAT and HESS. 



B. Secondary electrons and positrons 

Positrons in GCR are generally assumed to be of purely 
secondary origin, arising through the decay of pions and 
kaons produced in the interactions of GCR protons (and 
nuclei) with the interstellar medium (ISM) [2 . The neu- 
tral pions decay into 7-rays which then contribute to, 
if not dominantly constitute, the Galactic 7-ray back- 
ground. The charged pions on the other hand decay into 
neutrinos and muons, the latter subsequently decaying 
into electrons and positrons. Assuming that spatial and 



temporal variations in the GCR proton flux Jp and the 
ISM gas density njsM are small, the source density of 
these secondary background e~ and e+ is also homoge- 
neous, both in space and in time: 



^ISM 



msM c 



Eth: 



where dcrpp^e±+x/d£^ is the partial differential cross- 
section for production and jS c^l'i^ the velocity of the 
GCR. We can then integrate the Green's function for a 
single source over space and time to calculate 



J±{E) 



1 



47r \h{E)\ 



where x and £cr are defined by Eq. ( A2 ) of Ap pendix LAI 
^ is the diffusion length defined by Eq. (A3), and h ~ 



and 



0.1 kpc is the height of the Galactic disk. 

We calculate the flux of secondary background e~ 
e+ from the Solar-demodulated flux of GCR protons as 
derived from the BESS data and model the cross- 
sections according to Ref. [54]. The contribution from 
kaon decay is subdominant and is therefore neglected. 
The presence of He both in GCRs and in the ISM is 
taken into account by multiplying the proton contribu- 
tion by a factor of 1.2. Our results are in good agreement 
with Ref. [55 , taking into account the different diffusion 
model parameters and keeping in mind that convection 
and reacceleration have been neglected here. These fluxes 
are shown (dashed line) in the middle panel of Fig. [5] 
and are clearly a subdominant component which cannot 
acount for the deficit at high energies. 

Moreover, the positron flux is falling at all energies 
whereas the PAMELA data [l] clearly show a rise above 
a few GeV. One way this can be resolved is if there is 
a dip in the electron spectrum between ~ 10 and 100 
GeV. It has been suggested that Klein-Nishina correc- 
tions to the Thomson cross section for inverse Compton 
scattering [5^^ or inhomogeneities in the distribution of 
sources [31] can produce such a dip. However the for- 
mer would require a rather enhanced interstellar back- 
ground light (IBL) field [56], while the latter calcula- 
tion EI] assumes an incomplete source distribution (see 
Sec. [H]) and moreover adopts diffusion model parameters 
quite different from those derived from the measured nu- 
clear secondary-to-primary ratios [57] and the measured 
Galactic magnetic field and IBL [28l. 

The other, perhaps more straightforward possibility is 
to consider an additional component of GCR positrons 
with a harder source spectrum that results in a harder 
propagated spectrum and therefore leads to an increase 
in the positron fraction. 



C. Secondary accelerated electrons and positrons 

It has been suggested that acceleration of secondary 
produced through pp interactions inside the same sources 
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where GCR protons are accelerated, e.g. SNRs, can pro- 
duce a hard positron component [19 . We recapitulate 
here the essential formalism of diffusive shock accelera- 
tion [161 Hz] which yields the spectrum of the accelerated 
protons. This serves as the source term for calculating 
the spectrum of the secondary e^. 

The phase space density, /±, of secondary e~ and e+ 
produced by the primary GCR, both undergoing DSA, is 
described by the steady state transport equation: 



downstream upstream 



dx 



d_ 

dx 



ox 



Idii df± 



where q± is the source term determined by solving an 
analogous equation for the primary GCR protons. (Ide- 
ally we should solve the time-dependent equation, how- 
ever we do not know the time-dependence of the param- 
eters and can extract only their effective values from ob- 
servations. This ought to be a good approximation for 
calculating ratios of secondaries to primaries from a large 
number of sources which are in different stages of evo- 
lution.) We consider the usual setup in the rest-frame 
of the shock front (at x = 0) where ui {U2) and rii 
(712) denote the upstream (downstream) plasma velocity 
and density, respectively. The compression ratio of the 
shock r = Ui/u2 = n2/ni determines the spectral index, 
7 = 3r/(r — 1), of the GCR primaries in momentum space 
(note 7 = 2 + r). To recover 7 4.4 as determined from 
7-ray observations (see Table [l| we set r ^ 3.1. As noted 
earlier the theoretical expectation is however r = 4. 
For X 7^ 0, Eq. (12) reduces to an ordinary differential 



equation in x that is easily solved taking into account the 
spatial dependence of the source term 



^0 .^(p)e:cni/D(pp) ^ ^ 

for X > 0, 



(13) 



where the proton momentum should be distinguished 
from the (smaller) momentum p of the produced secon- 
daries, the two being related through the inelasticity of 
production: ^ :^ 1/20. Assuming D oc p (Bohm dif- 
fusion) in the SNR, the solution to the transport equa- 
tion p2| across the shock can then be written (see Ap- 
pendixIBI) : 



Tc/di 



ni \^ 



pCa=/dl 



for X < 0, 
for X > 0, 



(14) 

where di = D/ui is the effective size of the region where 
e~ and e+ participate in DSA (see Fig.j4|. 

The coefficients appearing in Eq. ( |14| satisfy an or- 
dinary differential equation dictated by continuity across 
the shock front (see Appendix [B|). This has the solution: 



fUv) 



Jo P' \P 



D{p')q±Ap') 



(15) 





FIG. 4: DSA setup in the rest frame of the shock front, ui 
{U2) and ni (712) denote upstream (downstream) plasma ve- 
locity and density, respectively. The right panel shows the 
solution of the transport equation for the primary GCRs. Par- 
ticles within a distance D/u of the shock front participate in 
the acceleration process. 



Assuming Feynman scaling for the pp interaction, i.e. 
pdcTpp/dp (X Ti± we can express the momentum depen- 
dence of the source term as 

/ X Cngas,l . /^.r / /x^CTpp^ei+X 



4:7rp' 
47rp^ 



dp 



Ncr{p) 



7-2 



(16) 



The maximum energy of protons is determined from the 
average maximum 7-ray energy £^max — 20 TeV (see Ta- 
ble |l| through the inelasticity of the pp ^ 7 + X process 
as - 20TeV/0.15 ^ 100 TeV [38]. 



We can easily interpret the solution (14) in terms of 
power laws in momentum. The second term downstream, 
{q2/u2)x, follows the spectrum of the primary GCRs 
((X p~^) and describes the production of secondary e~ 
and e+ that are then advected away from the shock front. 
However, secondaries that are produced within a distance 
~ D/u from the shock front are subject to DSA (see 
Eq. [13] and Fig. [4]). The fraction of secondaries that 
enters the acceleration process is thus given by the ra- 
tio of the relevant volumes, i.e. {D / ui) / {u2 ts-nr) ^ and 
the number density injected into the acceleration process 
is (1/^ + r'^)Dq±^i/u\. This rises with energy because 
of the momentum dependence of the diffusion coefficient 
{D{p) (X p) so the first term downstream in Eq. (14) gets 
harder: f±{p) (xp~^+-^. 

The injection spectrum R± is obtained by integrating 
the steady state solution over the volume of the SNR: 

i^±=4V47r/ dxx'^f±{x,p). (17) 

The resulting source spectrum, R±^ is thus the sum of 
two power laws. 



(18) 



7+2 


1 + 


( ' )\ 






VPcross / _ 



R± c^Rlp- 

where the "cross-over" momentum, Pcross, satisfies 

3 rulrsNR 



-^(Pcross) — 



4 7(l/e + r2) ■ 



(19) 
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As has been noted [19], this mechanism is most efficient 
for old SNRs where field amphfication by the shock wave 
is not very effective anymore. We therfore introduce 
a fudge factor Kb that parameterises the effect of the 
smaher field amplification on the otherwise Bohm-like 
diffusion coefficient in the SNR, 



D{E) = 3.3 X 10^^ Kb 



B 

'JIg 



E \ 2-1 



(20) 



The number of particles entering the acceleration 
process can of course not exceed the total number of 
secondaries produced inside the SNR. This effectively 
caps the growth of the term D{p')q^ i{p^)/u1 once 
{D/ui)/{u2rs^R) becomes larger than unity, a relation 
that defines a characteristic momentum scale Pbreak- We 
therefore substitute in Eq. (15), 



i^(Pbreak)g±,i(p) 



for p < Pbreak, 
for p > Pbreak. 



(21) 

The source spectrum R± thus returns to a p ^ depen- 
dence around p = Pbreak- At even higher energies the 
secondary spectrum cut s off a t the same £^cut as for pri- 
mary electrons (see sec. Ill A). 

Following Refs. [19, 23 , the parameters are chosen to 
be: ui = 0.5 x lO^cms"^, ngas,i = 2cm~^, B = 1 /iG. 
Choosing r = 3.1 to recover F = 2.4 the characteristic 
momenta Pcross and Pbreak turn out to be. 



(22) 
(23) 



What is still missing is the normalization of the injec- 
tion spectrum, i?^, in the sources which is proportional 
to the normalisation of the OCR protons, A^gcr, through 
Eq. (IT3|. Usually a factor K^p ^ 10'^ - 10"^ is intro- 



duced to normalize the electron component with respect 
to the protons; this depends on how particles are injected 
from the thermal background into the acceleration pro- 
cess and is not reliably calculable from first principles. 
We can get around this by assuming that the 7-rays de- 
tected from known SNRs by HESS are of hadronic origin, 
as is expected in this framework. Thus we can use the 
total luminosity of individual sources in 7-rays, 



(24) 



to determine the normalization of the proton component 
and therefore also the secondary injection rate g^. 

The compilation of 7 ray data on SNRs from HESS, 
see Table ^1 suggests an average value :^ 5.7 x 
10^^ s~^ TeV~^. We find then for the total spectrum 



(25) 



primary e 

















10^ 

Energy [GeV] 



secondary e + e' 



10^ I 




10^ 10' 
Energy [GeV] 




FIG. 5: Predicted spectra of electrons and positrons with data 
from Fermi LAT 4 (red circles) and HESS (blue squares 
& green triangles) . The diagonal arrows show the energy scale 
uncertainty. Top: Primary electrons after propagation to 
Earth. Middle: Secondary electrons and positrons from cos- 
mic ray interactions, created during propagation (dashed line) 
and created during acceleration in SNRs (full lines). Bottom: 
The sum of primary and secondary electrons and positrons. 



where 11+ (and analogously S^) are defined by Eq. (16), 
or explicitly 



^lOVAs.Tx 1033s-iTeV-V^''^ ' 

(26) 
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Diffusion Model 



Do 


10^^ cm^ s" 


1 


1 from GCR nuclear 


S 


0.6 




1 secondary-to-primary ratios 


L 


3 kpc 






b 


10-'^ GeV-' 




CMB, IBL and B energy densities 


Source Distribution 


^max 


1 X 10^ yr 


from ^nain ^ 3.3 GcV 


TSNR 


10^ yr 




from observations 


N 


3 X 10^ 




from number of observed SNRs 



Source Model 





1.8 X 10^°GeV-^ 


fit to e flux 


at 10 GeV 


r 


2.4 


average 7-ray spectral index 




20TeV 


typical 7-ray maximum energy 


Ecut 


20TeV 


DSA theory 




K 


7.4 X lO^^GeV-^ 


cf. Sec. 


IIIC 




Kb 


15 


only free paramet 


er (for fixed F) 



TABLE II: Summary of parameters used in the Monte Carlo 
simulation, for an injection spectral index F ^ 2.4. 



In the Monte Carlo code we have explicitly input the 
experimentally measured pp cross-section which gives a 
similar normalisation as the estimate presented above 
assuming Feynman scaling. The normalisation for sec- 
ondary electrons is computed similarly. 

The middle panel of Fig.[5]shows an example of the flux 
of secondary source e~ and e+ for 30 "histories" of SNRs 
in our Galaxy. Clearly this component can potentially 
match the high energy Fermi LAT and HESS data. 

We note that in our model, the contribution from sec- 
ondary electrons and positrons to the total flux is about 
twice as large as in ref.[19 where the primary injection 
spectrum was assumed to be ex motivated by DSA 
theory. However this is not consistent with 7-ray obser- 
vations of SNRs as seen from Tabled! 



IV. RESULTS 

The parameters used in the Monte Carlo are given in 
Table [Hj For an assumed injection spectral index F, the 
only free parameter is Across {cf- Eq. (Il9|) or, equivalently. 



the factor Kb {cf. Eq. (20)) which is determined by fitting 
the total flux of electrons and positrons to the Fermi 
LAT and HESS data (see Fig.[5|. Adopting F = 2.4, we 
find good agreement for Kb — 15, which corresponds to 
a cross-over of the primary and (accelerated) secondary 
components at Pcross — 28 GeV and a spectral break at 
Pbreak ^ 510 GcV {cf. Eqs. ^ and ([23|). 

We have calculated the x^~^ with respect to the com- 
bined Fermi LAT and HESS data for each configura- 
tion m of source distances and times, over 
all energy bins j. The three best "fits" are shown in 
Fig. [6] for different values of Kb and for F = 2.3, 2.4 
and 2.5 (see Table [Il|) . The corresponding predictions for 



the e+ fraction are shown in the bottom panels. These 
agree reasonably well with the data down to 6 GeV; we 
would not expect agreement at lower energies since we 
have neglected convection and reacceleration during in- 
terstellar propagation. In fact the PAMELA measure- 
ments of the e+ fraction are systematically lower than 
previous measurements, e.g. AMS-01 or HEAT, and it 
has been noted that this discrepancy can be resolved by 
considering charge-sign dependent Solar modulation with 
^+ = 438 MV for e+ and ^_ = 2MV for e' [52^ (rather 
than 0+ = ^_ = 600 MV). This however seems to be 
at odds with preliminary PAMELA data on the abso- 
lute electron flux [50 which does show substantial Solar 
modulation. Accordingly in Figj6]we have shown the pre- 
dicted e+ fraction for both cases; note that this does not 
affect our predictions for energies above 10 GeV. 

Thus our fits to both the PAMELA and the Fermi LAT 
spectra, including secondary e+ accelerated in SNRs, 
provides a consistent picture of current data on cosmic 
ray e~ and e+ between a few GeV and tens of TeV. Turn- 
ing the argument around, since a large fraction of the e~ 
and e+ observed in GCR above hundreds of GeV are re- 
quired to be secondaries in this model, there must be a 
large number of hadronic cosmic ray accelerators in our 
Galaxy, some of which should be quite nearby. 

An independent test of the model is provided by the 
usual 'messengers' of such hadronic acceleration environ- 
ments, namely 7-rays and neutrinos. Taking the known 
distribution of SNRs in the Galaxy (see Sec. [ll| we have 
calculated the column depth in SNRs in the Galactic disk 
as seen from Earth, 



POO 

X{4>')= / &r'r'g{r{r', <!>'), <l>{r',4>')), 
Jo 



(27) 



and show this in the top panel of Fig. [7| As expected, 
the column depth is largest towards the Galactic centre. 
However, the quantity that is more important for obser- 
vations is the brightness of sources. We have therefore 
weighted the integrand in Eq. (27) by 1/r^ and this flux 
weighted column depth is also shown in the top panel of 
Fig. [71 We note that although the maximum brightness 
is still expected around the Galactic Centre, the bright- 
ness in other directions is smaller by only ~ 30% because 
the sources in the closest spiral arms are then dominant 
(if they are actually there of course). 

This is illustrated in the bottom panel of Fig. [7] by 
an example distribution of SNRs from the Monte Carlo 
simulation, denoted by circles. The position of the cir- 
cle denotes the Galactic longitude and the radius is pro- 
portional to the brightness in units of the Crab Nebula, 
i.e. an integrated flux of (1.98 ± 0.08) x 10"^^ cm-^ s"^ 
above 1 TeV [58] . For a source of luminosity of = 
5.7 X 10^^TeV"^s-^ (see Table |l| at distance d, the in- 
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FIG. 6: The three best fits (out of 30 source "histories") to the total spectrum of electrons and positrons measured by Fermi 
LAT [4] (red circles) and HESS [51 16] (blue squares & green triangles), and the corresponding prediction for the positron fraction 
for different values of F and Xb, for both charge-sign independent (full line) and charge-sign dependent (dashed line) Solar 
modulation (see text for details). The PAMELA data is shown for comparison (open red circles). 
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FIG. 7: Top: The column depth and flux weighted column 
depth of the SNR density in the Galactic plane. Bottom: 
Example of a distribution of SNRs in 7-rays/neutrinos from 
the Monte Carlo simulation. The position of a circle denotes 
the Galactic longitude of the source and the radius is pro- 
portional to the brightness in units of the Crab nebula. One 
source whose circle exceeds the vertical scale is ^ 500 pc from 
Earth and has a total integrated flux above 1 TeV of 6 times 
the Crab Nebula. 



tegrated flux above 1 TeV is, 
1 



-2 



cm ^ s 



(28) 



i.e. about 40 % of the Crab Nebula flux at (i = 2 kpc. It 
is seen that although most of the sources are clustered to- 
wards the Galactic centre, there are several bright sources 
at large longitudes as well. We find typically ~ 3 sources 
brighter than the Crab (or ~ 7 brighter than 50 % Crab) . 

The adopted distribution of SNRs (Sec. [ll| and the 
average luminosity per source determined from a compi- 
lation of known sources (Table [l|) thus leads to the pre- 
diction of several nearby SNRs with fluxes of the order of 
the Crab Nebula. Note, however, that close sources could 
be rather extended and thus have escaped detection by 
HESS in one of its surveys of the Milky Way [44", ^ [61] . 
For example, a diameter of ~ 50 pc which is a typical 
value for a very old SNR, corresponds to 1.5° at 2 kpc. 

Extended 7-ray luminous SNRs can however be de- 
tected by MILAGRO [59^ with its larger field of view. A 
survey in Galactic longitude / G [30°, 220°] and latitude 
h G [—10°, 10°] has revealed 6 new sources at a median 



energy of 20 TeV, several of which are spatially extended. 
The flux from a SNR of the above luminosity at <i = 2 kpc 
is Qy{AiTdP) ^ 1.2 X 10-11 TeV"^cm-2s-i at 1 TeV. 
Scaled with a spectral index of 2.4 to 20 TeV, this gives 
Q^/(47rd2) 20-2 '^ c:^ 9.0 x 10-^^ TeV"i cm-^ s-^ which is 
in the range of the unidentified MILAGRO sources [59] . 
We note that the MILAGRO source MGRO J1908+06 
was recently confirmed by HESS [62 , though with a 
smaller angular extent of ~ 0.7°. However, correlating 
unidentified MILAGRO sources with the FERMI BSL 
list [63l seems to favour associations with pulsars, al- 
though several new unidentified extended sources have 
also been found. 

Hadronic sources of cosmic rays should also be visi- 
ble by their neutrino emission. On general grounds, the 
neutrino luminosity (from tt^ decay) can be directly re- 
lated to the 7-ray luminosity (from tt^ decay) and should 
be of the same order of magnitude since pp interac- 
tions produce 7r+, tt^ and tt" in roughly equal numbers. 
Each of the three neutrinos produced in the decay chains 
7r+ ji^v^ e+i^e^/iZ^/x and tt" fi'i^f^ ^ e~Uei^^iy^ 
carries about half of the energy of each photon produced 
in the decay tt^ 77. Hence, the ratio of neutrinos to 
photons produced on average is ~ 3 : 1 and the total 
neutrino luminosity is 



r^o / 

^ ^ TeVy 



(29) 



Presently the largest cosmic neutrino detector is the Ice- 
Cube observatory [65] under construction at the South 
Pole. IceCube observes high energy neutrinos via their 
interactions with nucleons in the vicinity of the detec- 
tor and subsequent Cerenkov light emission of energetic 
charged particles in the transparent glacial ice. The most 
important signal for neutrino astronomy is the Cerenkov 
radiation by muons produced via charged current in- 
teractions of muon neutrinos. Since the muon inherits 
the large boost of the initial neutrino the point source 
resolution is ~ 1°. The large background signal of at- 
mospheric muons is efficiently reduced for upward-going 
muons, i.e. neutrino sources which are somewhat below 
the horizon. Hence, IceCube is mainly sensitive to neu- 
trino point sources in the northern sky, which excludes 
SNRs in the direction of the Galactic centre. 

Neutrino emission associated with galactic TeV 7-ray 
sources has been investigated by many authors [66l [671 
l68] l69] [70] [71] including also the HESS sources used 
in our analysis. In particular, Ref. [67 investigates 
the prospects of neutrino detection for the SNRs HESS 
J0852.0-463, J1713-381, J1804-216, J1834-087 (see Ta- 
ble [l| in the proposed KM3NeT detector in the Mediter- 
ranean which will see the Galactic centre region. The 
muon neutrino rate is expected to be a few events per 
year for such sources. 

Due to flavour oscillations of neutrinos with large mix- 
ing angles, the initial flavour composition Qj^^ : Qjy^ : 
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Qjy^ :^ 1 : 2 : from pion decay is expected to be- 
come ~ 1 : 1 : 1 at Earth. The TeV muon neutrino 
point flux from a hadronic 7-ray source located at a dis- 
tance d and with a power-law index T 2.4 is thus 
F^^{> 1 TeV) 2^-^F^{> 1 TeV), hence 



F^^{> 1 TeV) 3.2 x 10" 



2 kpc 



cm ^s ^. 



(30) 



This should be compared to the results of searches for 
neutrino point sources in the northern sky, in particu- 
lar the close-by SNR Cassiopeia A (see Table [ij), using 
data taken with AMANDA-II (the predecessor of Ice- 
Cube) during 2000-2006 [72 and, more recently, with 
the first 22 strings of IceCube during 2007-08 [73 . The 
average 90% C.L. upper limit on the integrated z^^ flux 
in the energy range 3 TeV to 3 PeV is [73] 



< 4.7 X 10" 



-12 -2 - 

cm s 



(31) 



i.e. well above the flux of ~ 7 x 10~ cm~^ s~^ expected 
from a SNR at 2 kpc, assuming T = 2.4. 

The full 80 string configuration of IceCube thus has ex- 
cellent prospects to identify these SNRs. A point source 
in the northern sky with an E~'^ muon neutrino flux. 



7.2 X 10" 



(32) 



in the TeV-PeV range can be detected with a 5<j sig- 
nificance after three years of observation. This does de- 
pend somewhat on the spectral index and energy cut-off, 
since the signal (after "level 2 cuts") peaks at an en- 
ergy of ~ 10 TeV [65]. As mentioned previously, our 
analysis predicts on average ^ 3 nearby 7-ray sources 
stronger than Crab with corresponding muon neutrino 
fluxes larger than ~ 7 x 10~^^ cm~^ s~^. Note that al- 
though the Galactic centre is not in the field of view of 
IceCube, SNRs following the spiral arm structure of the 
Galaxy are expected to be detected also in the Galactic 
anti-centre direction, as seen in the example distribution 
shown in the bottom panel of Fig. [7| 



V. SUMMARY 

Supernova remnants have long been suspected to be 
the sources of Galactic cosmic rays. We have discussed 
a recent proposal [19 that proton-proton interactions in 
the shocks of SNRs followed by the diffusive shock accel- 
eration of the secondary positrons produced can flatten 
the spectrum of the secondaries relative to that of the 
primaries. These hard spectra may be the origin of the 
recently observed cosmic ray "excesses" — both the e+ 
fraction observed by PAMELA [T] and the e~ + e+ flux 
measured by Fermi LAT [4] and HESS [5l|6]. 

We have investigated how 7-ray emission of SNRs - as- 
sumed to be of the same hadronic origin as the positrons 



- together with cosmic ray data, constrain the acceler- 
ation of positrons. We have accounted for the spatial 
and temporal discreteness of SNRs via a Monte Carlo 
exercise, drawing samples from a realistic galactic distri- 
bution with the observed SN rate. For the diffusion pa- 
rameters we have adopted standard values derived from 
cosmic ray nuclear- to-primary ratios, as well as the en- 
ergy densities of galactic radiation and magnetic fields. 

We have compiled a list of all 7-ray emitting SNRs 
observed by HESS and determined the mean value of 
the flux, which fixes the hadronic interaction rate in the 
SNR. Low energy data from PAMELA on the absolute 
e~ flux was used to normalize the primary flux of e~ . 
The contribution from accelerated e+ was then found by 
fitting the e~ + e+ flux to Fermi LAT and HESS data, 
adjusting the (only) free parameter Kb which determines 
the diffusion rate near SNR shocks. 

The spectra of e+ and e~ thus derived agrees well 
with the e+ fraction observed by PAMELA in the range 
5 — 100 GeV. The apparent deficit at lower energies can 
be attributed to convection and diffusive reacceleration 
of primary electrons that become important at these en- 
ergies and were neglected in our analysis. The flux of e+ 
and e~ becomes dominated by the accelerated secondary 
component at high energies; the corresponding e+ frac- 
tion levels out at ~ 0.4, reflecting the relative multiplicity 
of e+ and e~ produced by p-p interactions. 

We note that although the value of ~ 10—20 found here 
for Kb differs from the value of ~ 40 we had found earlier 
from our analysis of the titanium-to-iron (Ti/Fe) ratio in 
cosmic rays [24], the latter determination is subject to 
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FIG. 8: The B/C ratio in cosmic rays along with model pre- 
dictions (after Ref. [24]) — the leaky box model with produc- 
tion of secondaries during propagation only (dot-dashed line), 
and including production and acceleration of secondaries in 
a nearby SNR (solid lines) for values of the diffusion coeffi- 
cient near the shock wave which best fit the spectrum (see 
Fig|6|). The dashed line corresponds to the value of the dif- 
fusion co-efficient required to fit the ATIC-2 data on Ti/Fe 
(from Ref. [24]), along with the la and 2a error bands. The 
data points are from HEAO-3-C2 (circles) i;^, ATIC-2 (tri- 
angles) and CREAM (squares) [76]. 
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large uncertainties due to poor experimental statistics at 
high energies. This affects our previous prediction [24 for 
the boron-to-carbon (B/C) ratio and we show in Fig. [s] 
the corresponding la and 2<j error bands, along with the 
prediction taking Kb ~ 10 — 20 as indicated by our new 
fits to the spectrum. The predicted upturn at high 
energies wih soon be tested by PAMELA and AMS-02. 
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To be consistent with our overall framework the 7-rays 
observed from SNRs have been assumed to be of hadronic 
origin. The known spatial distribution of SNRs then im- 
plies (on average) several nearby sources with a 7-ray 
flux comparable to the Crab. We have speculated that 
some unidentified MILAGRO sources [59] might corre- 
spond to such old SNRs. Moreover, the same hadronic 
processes in SNRs will inevitably produce high energy 
neutrinos which can be detected in cubic-km telescopes 
such as IceCube [65 . The neutrino luminosity can be 
directly related to the 7-rays and is not connected to 
the hypothetical acceleration of e+ and e~ in the sources 
as in our present model. Nevertheless, similarly to the 
previous argument, we expect on average a few nearby 
sources, some of which may may also lie within the field 
of view of IceCube and can thus be detected with high 
statistical significance after three years of data taking. 



While our calculational framework is based on first- 
order Fermi acceleration by SNR shock waves, we have 
noted that in detail the observations do not fit the the- 
oretical expectations, e.g., the shock compression ratio 
inferred from the observed 7-ray spectrum (~ E~'^-^) 
is 3.1 rather than 4 as is expected for a strong shock 
[16]. Going beyond the test particle approximation, the 
generic expectation in such a process is for particle spec- 
tra which are much flatter than those observed (~ E~'^-^ 
and slightly concave)^ when the back reaction of the cos- 
mic rays on the shock is taken into account [17 . By 
contrast, the observed radio spectrum of Cassiopeia A is 
slightly convex and this, as well as the morphology and 
time evolution of radio emission from such young SNRs, 
can be well explained in terms of second- order Fermi ac- 
celeration by plasma turbulence behind the shock wave 
[77] . Moreover the observed spatial correlation between 
the 7-ray emission and the hard X-ray emission from 
some SNRs argues for a leptonic rather than hadronic 
origin and further observations are necessary to resolve 
this issue [78 . It has been argued that cosmic ray protons 
and nuclei may well have different sources (e.g. "super- 
bubbles" formed by multiple supernovae) than the cos- 
mic ray electrons [79^. The additional predictions made 
in this paper concerning the visibility of hadronic acceler- 
ators in 7-rays and neutrinos, tied to the expectations for 
the fluxes of the accelerated secondary positrons in cos- 
mic rays, will hopefully enable further consistency tests 
of the SNR origin hypothesis for galactic cosmic rays. 



APPENDIX A: GREEN'S FUNCTION OF 
DIFFUSION EQUATION 

We solve Eq. ([T]) for a halo of extent ±L in z direction, 
neglecting the boundaries in the radial direction. The 
Green's function for the flux of electrons from a source 
at r that went off a time t ago with a spectrum Q{E), is: 
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(7r£2)! 
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where 



V'^* n=-oo 

and the diffusion length I is defined as 
^2 - 4 / AE' > 
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with £cr = 4L/7r. If we neglect the spatial extent of the 
disk and set 2: = 0, the function x{^^^) is approximately: 



x(o,0 



-e 

1 



for / > 
for / <C 



(A4) 



In practice both limits can be connected at / :^ 0.66 
such that the approximated x{^^^) has a relative error 
of at most 0.5%. We motivate the choice of the param- 
eters of our diffusion model from an analysis of nuclear 
secondary-to-primary ratios [STj: I^o = lO^^cm^s"^, S = 
0.6, L = 3 kpc, and from the Galactic magnetic field and 
interstellar radiation fields [28,: 60 = 10"^^ GeV"^ s'^. 



APPENDIX B: BOUNDARY CONDITIONS OF 
DSA 



The solution (14) of Eq. 12 satisfies 

lim /± = , lim — — = and 

— — ox 



hm f± 



< 00 , 
(Bl) 
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Continuity at the shock front a; = requires: 



D 



df± 

dx 



D 



x=0- 



df± 

dx 



x=0+ 



1, , dfl 



yielding the differential equation, 



dp 



This is readily integrated with boundary condition 
f±{0) = and yields the p dependence in Eq. (15). 



(B2) 
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